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1 Introduction 

In f4 , following up work of Hitchin [9], the author found it useful to express 
Nahm's equations, for a matrix group, in terms of the motion of a particle in 
a Riemannian symmetric space, subject to a potential field. This point of view 
lead readily to an elementary existence theorem for solutions of Nahm's equa- 
tion, corresponding to particle paths with prescribed end points. The original 
motivation for this article is the question of formulating an analogous theory for 
the Nahm equations associated to the infinite-dimensional Lie group of area- 
preserving diffeomorphisms of a surface-in the spirit of We will see that 
this can be done, and that a form of the appropriate existence theorem holds — 
essentially a special case of a result of Chen. However the main focus of the 
article is not on existence proofs but on the various formulations of the problem, 
and connections between them. In these developments, one finds that the nat- 
ural context is rather more general than the original question, so we will start 
out of a different tack, and return to Nahm's equations in Section 5. 

Consider the following set-up in Euclidean space R^, in which we take co- 
ordinates {xi,X2,z) — thinking of z as the vertical direction. (We will use the 
notation = ^iT'§^ — ^z-) Suppose we have a strictly positive function 
H{xi,X2)- This defines a domain 

^Ih ^ {{xi,X2,z) : < z < H{xi,X2)}, 

whose boundary has two components {z = 0} and {z ~ H}. We consider the 
Dirichlet problem for the standard Laplacian: to find a harmonic function 9 on 
ri// with 9 — on {z = 0} and 6* = 1 on {z = H}. To set up this problem 
precisely, let us assume that the data H is Z^-periodic on R^, and seek a Z^- 
periodic solution 9. Now we have a unique solution 9 to our Dirichlet problem. 
Consider the flux of the gradient of 9 through the boundary {z — H}. This 
defines another function p on R^. To be precise, if lh is the obvious map from 
R^ to the boundary {z = H} then the fiux is defined by 

i*u{*d9) = pdxidx2- (1) 
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Explicitly 



p = d,e- {di0diH + d2ed2H), 



with the right hand side evaluated at (xi, 2:2, i/(a;2, 2:2)). By the maximum 
principle, p is a positive function, since the normal derivative of 9 is positive in 
the positive z direction on {z = H}. We consider the following free boundary 
problem: given a positive periodic function p does it arise from some periodic 
H, and is H unique? 

One can gain some physical intuition for this question by supposing that 
that the lower half-space {z < 0} represents a body with an infinite specific 
heat capacity fixed at temperature and corresponds to a layer of ice 
covering this body. We choose units so that the melting temperature of the ice 
is 1. Sunlight shines vertically downwards onto the upper surface {z = H} of 
the ice, but with a variable intensity so that heat is transmitted to the surface 
according to the density function p. We suppose that the surface of the ice 
is sprinkled by rain, which will instantly freeze if the surface temperature of 
the ice is less than 1. We also suppose that a wind blows across the surface, 
instantly removing any surface water. Then we see that the solution to our 
free boundary problem represents a static physical state, in which the upper 
surface of the ice is just at freezing point, the lower surface is at the imposed 
sub-freezing temperature and the heat generated by the given sunlight flows 
through the ice without changing the temperature. Physical intuition suggests 
that there should indeed be a unique solution. 

We can express the free-boundary problem considered above as a special 
case of another question. Suppose now that we have a pair of periodic functions 
Hq, Hi on with Hq < Hi. Then we have a domain QhoMi = {Ho{xi,X2) < 
z < Hi{xi,X2)}, with two boundary components. Let 9 be the harmonic func- 
tion in this domain equal to 0, 1 on {z = Ho},{z = Hi} respectively. Then 
we obtain a pair of flux- functions po,Pi as before. By Gauss' Theorem, these 
satisfy a constraint 



since [0, 1]^ is a fundamental domain for the Z^-action. Obviously, if we replace 
Hq, Hi by Hq + c. Hi + c for any constant c we get the same pair po) Pi- We ask: 
given Pq,Pi satisfying the integral constraint (2), is there a corresponding pair 
{H{),Hi), and if so is the solution unique up to the addition of a constant? A 
positive answer to this question implies a positive answer to the previous one, by 
a simple reflection argument. (Given p, as in the first problem, take po = Pi = 
p/2. Then uniqueness implies that the solution has reflection symmetry about 
9 = 1/2 and we get a solution to the flrst problem by changing 9 io 29 — 1.) 

Of course we can also imagine a physical problem corresponding to this 
second question: for example a layer of ice in the region ^Ha,Hi- We can 
now vary the problem by supposing that in place of ice we have a horizontally 
stratified material in which heat can only flow in the horizontal directions. Thus 
the steady-state condition, for a temperature distribution 6{xi,X2-iZ) is 





{dl + dl)9 = Q. 



(3) 
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We define fliix- functions po, Pi by pulling back the 2-forni 

di6 dx2dz — 826 dxidz, 

and the integral constraint (2) still holds. So we ask: given po,Pi satisfying 
(2), is there a pair Hq, Hi and a function 9 on Q,Ho,Hi , equal to 0, 1 on the two 
boundary components, which has these fluxes, and is the solution essentially 
unique? (In this case one has to relax the condition on the domain to Hq < Hi.) 

It is natural to extend these questions to a general compact oriented Rieman- 
nian manifold X (which would be the flat torus R^/Z^ in the discussion above). 
Write dp for the Riemannian volume form on X. We fix a real parameter e > 
and define a map *^ from T*{X x R) to A"T*(X x R) by 

*edz = edp , *ea = {*xct)dz, 

for a e T*X, where *x is the usual Hodge *-operator on X. Then for a function 
^ on a domain in X x R 

d*^d6= {A^9)dzdp, 

where 

A,^ = (-6^2 + Ax)e, 

with Ax the standard Laplace operator on X. (We use the sign convention 

that Ax is a positive operator, so when e = 1 our A^ is the standard Laplace 
operator on X x R.) If 9 is defined on a domain Qho,Hi, as above, we define 
the flux Pi on the boundary {z = Hi} by pulling back *fd9 just as before. We 
consider a pair of functions po,Pi > with 

Po dp = pi dp = dp 
Jx Jx Jx 

and we ask 

Question 1 Is there a pair Hq < Hi and a funetion 9 on the set ^IhoMi ^ 
X X H with 9 = 0,1 on the hypersurfaces {z = Hq}, {z = Hi}, with fluxes Pi 
and with A^9 = 0? If so, is the solution essentially unique? 

For any e > the equation A^^ = can be transformed into the standard 
Laplace equation on the product, by rescaling the z variable. When e = the 
equation has a very different character: it is not elliptic and we obviously do 
not have automatic interior regularity with respect to z. 

2 An infinite- dimensional Riemannian manifold 

We now start in a different direction. Given our compact Riemannian manifold 
X we let H be the set of functions ^ on X such that 1 — Ax(t> > 0. We make 
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H, into a Riemannian manifold, defining the norm of a tangent vector 5(j) at a 
point <f> by 

\2 _ (s:±\2 



\mi^ / {HY {i-^x4>)diji. 



X 



Tlius a path <j)[t) in 7i, parametrised by i e [0, 1] say, is simply a function on 
X X [0, 1] and the "energy" of the path is 

When X is 2-dimensional and orientable, this definition coincides with the met- 
ric on the space of "Kahler potentials" discussed by Mabuchi llj, Semmes [12] 
and the author [6] . The general context in those references is a compact Kahler 
manifold: here we are considering a different extension of the 2-dimensional 
case, and we will see that some new features emerge. The account below follows 
the approach in [5] closely. 

It is straightforward to find the Euler-Lagrange equations associated to the 
energy (3). These are 

These equations define the geodesies in Ti. We can read off the Levi-Civita 
connection of the metric from this geodesic equation, as follows. Let (j){t) be 
any path in H. and ^{t) be another function on X x [0, 1], which we regard as 
a vector field along the path Then the covariant derivative of (j) along the 
path is given by 

A0 = ^ + (m,Vx^), (5) 

where 

and ( , ) is the Riemannian inner product on tangent vectors to X. (We write 
Vx, or sometimes just V, for the gradient operator on X, so Wt is a vector 
field on X .) This has an important consequence for the holonomy group of the 
manifold H. Observe that the tangent space to 7i at a point is the space of 
functions on X endowed with the standard inner product associated to the 
measure 

= (1 - A0)d/xo. 

So, in a general way, the parallel transport along a path from (j)Q to 4>i should 
be an isometry from (X, d/i^o ) to {X, dficjii ) . (Here we are ignoring the 
distinction between, for example, smooth functions and functions.) What 
we see from equation (5) is that this isometry is induced by a diffeomorphism 
f : X ^ X with 

f*{dn^i) = dn^o- (6) 
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The difFeomorphism is obtained by integrating the time-dependendent vector 
field Wt and equation (6) foUows from the identity 

CwtdjJ,^ = d * (- — \r—d4> *x djjL^) = = 

i — Ax at 

(Here C denotes the Lie derivative on X.) We conclude that the holonomy 
group of H is contained in the group Q of volume-preserving diffeomorphisms of 
{X, dfio), regarded as a subgroup of the orthogonal group of L'^{X, dfio). (This 
can also be expressed by saying that there is an obvious principal ^-bundle over 
H with the tangent bundle as an associated vector bundle, and the Levi-Civita 
connection is induced by a connection on this C/-bimdle.) 

We now move on to discuss the curvature tensor of Ti. Let be a point of 
H and let a, /3 be tangent vectors to H at (p — so a and /3 are just functions on 
X. The curvature Ra.p should be a linear map from tangent vectors to tangent 
vectors: that is from functions on X to functions on X. The discussion of the 
holonomy above tells us that this map must have the form 

Ra.A^) = i''c,,(3,^^), (7) 

for some vector field iya,i3 on X, determined by </>, a, (3. Moreover we know that 
we must have 

Cu^^ {diJ.4,) = 0. 

To identify this vector field we introduce some notation. For vector fields v, w 
on X we write v x w for the exterior product: a section of the bundle A^TX. 
We define a differential operator 

curl : T{A^TX) T{TX), 

to be the composite of the standard identification: 

A'^TX ^ A"-2r*x, 

(using the Riemannian volume form dfj.), the exterior derivative 

d : T{A"-'^T*X A"-'^T*X, 
and the standard identification 

A"''^T*X = TX, 
(using the volume form dfj, again). Then we have 
Theorem 1 The curvature ofH is given by (7) and the vector field 

v^S = T^curl (^^Va x V/3). 

Corollary 1 The manifold H has non-positive sectional curvature. 
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The sectional curvature corresponding to a pair of tangent vectors a, (3 at a 
point (j) is 

In our case this is 



Ka,0= / (z^a,/3,Va)/3(l - Ax</')dAt. 
Jx 



Unwinding the algebraic identifications we used above, the integrand can be 
written in terms of differential forms as 

. \ J '^^df- i—^ {da A dp)] (3{1-Ax<t>). 



So 



Ko^.fi = I da Ad [y^Ta^ * ^'^'^ ^ '^^^ I ^^ 



Applying Stokes' Theorem this is 

Ka,f3 = - I - — ^^—-daAdl3A*{daAdl3) = - ( - — ^——-\daAdl3\^ djjL < 0. 

In the proof of Theorem 1 wo will make use of two identities. For any pair 
of vector fields v, w and function / 

curl {v X w) = [v, w] + (div v)'w — (div iu)v (8) 

curl {f{v X w)) = /curl {v x w) + {v, V f)w - {w, Wf)v. (9) 

Wc leave the verification as an exercise. (Considering geodesic coordinates we 
see that it suffices to treat the case of Euclidean space. Our notation has be 
chosen to agree with standard notation in the case of vector fields in R^.) 

To calculate the curvature wc consider a 2-parameter family (/)(s, t) in ?i, with 
a corresponding vector field tp{s, t) along the family. Then we will compute the 
commutator {DgDt — DtDs)'4'{s,t). Evaluating at </> = 0(0,0) this is Ra^p^ip) 
where ip = ip{0, 0), a = ds(j>, /3 = dtcj). 

Now we write 

Ds = j- + Ws , Dt = - + Wt, 
OS ot 

where the vector fields Ws,Wt are regarded as operators on the functions on X. 
So DgDt — DtDg is the operator given by the vector field 

and y is exactly the vector field ya,i3 we need to identiiy. Recall that 

Ws = T— , Wt - 



1 - ' ' 1 - A0' 
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So 

dW, -I d'^<p\ 1 



dt \dsdtj (l-A(/))2 

Evaluating at s = t = where ds(j) — a,dt(j) — (3 we have 

dWs dWt 1 



dt ds (1 - 



(AaV/3 - A/3Va) . 



Write g for the function (1 — A(p) ^. Combining with the Lie bracket term we 
obtain 

[gVa,gV(3] + g^{AaV/3 - A/3Va). 
Now applying (8) we have 

[gVa,gVf3] = curl (g^Va x V/3) + div {gVa)VP - div (gV/3)Va. 

Applying (9) we have 

curl (g^Va x V/3) 5 curl (gVa x V/3) + g{{\/g, Va)V/3 - (Vg, V^)Va). 

Since 

div igVa) = gAa - (Vg, Va) , div (.gV/3) = 5A/3 - (V5, V/3) 
we see that 

z^a,/3 = gcurl (^Va x V/3), 

as required. 

In the case when X has dimension 2 — as discussed in [6j, [11], [12] — the 
space Ti, is formally a symmetric space. This is not true in general, since the 
curvature tensor is not preserved by the action of the group Q. 

We define a functional on H by 

1/(0) = / ^d^l. 

J X 

This function is convex along geodesies in 7Y, since the geodesic equation implies 
> 0. Now introduce a real parameter e > and consider the functional on 
paths in Ji: 

j\\^\l + ^V{^) dt, (10) 

corresponding to the motion of a particle in the potential —eV. The Euler- 
Lagrange equations are 

■■ \Vx^? + e 
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3 Three equivalent problems 

In this section wc will explain that there are three equivalent formulations of 
the same PDE problem associated to a compact Riemannian manifold X. We 
have essentially encountered two of these already. 

• The "6* equation" . 

This is the problem we set up in Section 1. We are given positive functions 
Po,pi on X, with 

[ Pidfi= [ dfi. (12) 
Jx Jx 

We seek a domain flHo,Hi C X x R defined hy Hq,H\ : X — »• R and a 
function 9 on ^Ha,H-i , equal to 0, 1 on the two boundary components, with 
fluxes Pq,Pi and satisfying the equation 

= 0. 



• The equation" 

Here we are given ^o^i on X, with 1 — A0j > 0. We seek a function 
$ on X X [0, 1], equal to (?io, on the two boundary components, with 
1 — A$ > for all t and satisfying the nonlinear equation 

g^(l-A.v*)-|v(^)P = . (13) 

As we have explained in Section 2, this is the same as finding a path in the 
space 7i, with end points 4>o,4'ij corresponding to the motion of a particle 
in the potential —eV. 

Now we introduce the third problem. 

• The "C/ equation" 

We are given positive functions (l)o,(i>i, with 1 — A^j > 0, as above. Define 
a function L on X x R by 

L{x, z) — max(0o(a;) — 4>i{x) + z, 0). 

We seek a function U (x, z) on X x R with U > L everywhere and 
satisfying the equation 

AM = (1 - A<Ao) (14) 
on the open set O where U > L. 
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The equivalence of these three problems (assuming suitable regularity for 
the solutions in each case) arises from elementary, but not completely obvious, 
transformations. We describe these now. 

• ^-equation ^-equation 



Suppose we have a solution ^ on a domain Qho,Hi- Then dzO = ^ 
is positive on the boundary components of Qho,Hi- The function dzO 
satisfies the equation Ae{dzO) = and it follows from this that dz9 is 
positive throughout the domain. This implies that, for any t G [0,1], 
the set 0~^{t) is the graph of a smooth function ht on X. By definition 
ho = Ho and hi = Hi. We also write this function as h{t,x) where 
convenient. For each fixed t we can define a function pt on X by the flux 
of *f^d9, just as before. 

We claim that 

We show this by direct calculation (there are more conceptual, geometric 
arguments). For simplicity we treat the case when the metric on X is 
locally Euclidean, so Ax = —^df where dt = ^ , for local coordinates 
Xi. The identity 

e{x,ht{x))=t 

implies that 
and 

Now 



diO + d^e dih = (16) 

d^e dth = 1. (17) 

^xht = - Y^i^i + {dih)dz)dih, 



and this is 



which is 



On the other hand the flux pt is given by pulling back the differential form 
*td9 on the product by the map x i— > {x, ht{x)) and this gives 

Pt = ed.e + ^^Y.^dief. 



So 
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This is 

So we see that dtPt = —Axht, since ed^dzO = — ^ • 9,9,^. 

Now the normahsation (13) implies that there is a function (po on X such 

that po = 1 — Ax4>o- For f > we define 4>t by 

t 

hrdr. 

We can also regard this family of functions as a single function $ on 
X X [0, 1], Then (15) implies that pt = l — Axcjjt for each t. We have 



and 



So 



Now since 



i 

^(l-A.a>) = . + ^( 

= d,ht = -^diO 



we can write the above as 

^(l-A.$) = e+|v4$|^ 

as required. 
• ^-equation =J> [/-equation. 



Here we suppose we have a solution ^{x,t) of the $ equation and we write 

<I>(.T,0) = 00, $(.x, 1) = 4>i. Wc essentially take the Lcgcndrc transform in 
the i-variable. The discussion is slightly more complicated when e = 0, 
so for simplicity we treat the case when e > and 9j $ is strictly posi- 
tive. Write Hi{x), H2{x) for the derivatives dt^ evaluated at (x, 0), (.t, 1) 
respectively, so Hq < Hi. We calculate first in the open set Qho,Hi- For 
each fixed x G X and each z in the interval {Ho(x) , Hi{x)) there is a 
t = t{x, z) such that z = dt^. We set 

U{x, z) = ^{x, 0) - $(a;, t) + zt. 

This defines a function U in rLHo,Hx- We define U outside this set by 
setting U{x,z) = Q ii z < Hq{x) and U{x,z) = L{x,z) = (j)o — (f)i — z if 
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z > Hi{x). It follows from the definitions that U is C^, that U > L and 
that the set where U > L is exactly ^Iho,Hi- We calculate on this set. 
Then d^U = t and 

dlU={d^^)-\ (18) 
Differentiating with respect to the parameters Xi we have 

diU = di<j>o - dii^, 

and 

Z>0 



dxi dtdxi 

Differentiating the identity z = 9t$ gives 

0= 7TT7^ + 



dtdxi dt'^ dxi ' 
so we can write 

Summing over i and using the formula (18) for we obtain 

1 

and so 



edlU - AxU =g2^{e+ |Vxt$r) - Ax<^o + 1, 



A,U=1- Ax(t>o. 



• [/-equation ^-equation 

Now suppose we have a solution U of A^U = po in a domain i^Ho,Hi, 
satisfying the appropriate boundary conditions, where po = 1 — Ax(t>o- 
We set 

oz 

Then A^0 = and = 0, 1 on the two boundary components. We have 
to check that the fluxes of *^dO on the boundary components are pi = 
1 — Ax(j)i- Consider flrst the boundary component where z = Hq. The 
flux is 

Now we have identities 

{diF){x,Ho{x)) = , {d,F){x,Ho{x)) = 0. 
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Differentiating the first of these with repect to Xi we get 



dfF + diHodid.F = 0, 
on the boundary. Differentiating the second gives 

d,d,F + d.Hod^F = 
on the boundary. Combining these we have 

{d^diFf = {dlF){d^F). 

Hence the flux is 

edlF + Y,d^F = po. 

i 

The argument for the other boundary component {z = -ffi(x)} is similar. 

4 Existence results and discussion 

We have set up a class of PDE problems associated to any compact Riemannian 
manifold, and seen that these have three equivalent formulations. In this section 
we will make some remarks about existence results, and comparison with the 
free-boundary literature. This discussion is unfortunately rather incomplete, 
mainly due to the authors limited grasp of the background. 

4.1 Monge- Ampere and the results of Chen 

For a function $ on X x (0, 1) write for the nonlinear differential operator 

g($)=a,2$(l-Ax$)-|Vx^$p. 

So our "^-equation" is g($) = e. When X has dimension 1 — a circle with local 
coordinate x — we can write Ax = —d^. and the equation is the real Monge- 
Ampere operator 

When X has dimension 2 the operator can be expressed as a complex Monge- 
Ampere operator. That is, we regard X as a Riemann siirfacc and identify the 
Laplace operator on X with idd. We take the product with a circle, with co- 
ordinate a, and let t = t + iahe a, complex coordinate on the Riemannn surface 

X (0, 1). Then, in differential form notation, our nonlinear operator is given 

by 

{ujQ + idd^f = q{^)uJodTdrr, 

where loq is the Riemannian area form of X lifted to X x x (0,1). Our 
Dirichlet problem becomes a Dirichlet problem for S^-invariant solutions of this 
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complex Monge- Ampere equation on X x x (0, 1). This was studied by Chen 
[2] and it follows from his results that, for any e > there is a unique solution to 
our problem, and hence an affirmative answer to Question 1 in this case. (Chen 
does not state this result explicitly, but it follows from the continuity method 
developed in [5], Section 3, that for any strictly positive smooth function ly on 
X X [0, 1] there is a solution of the equation q{^) — v with prescribed boundary 
values (^0, 01-) 

It seems quite likely that the techniques used by Chen can be extended to 
the higher dimensional case. The foundation for this should be provided by a 
convexity property of the nonlinear operator which we will now derive. Let A 
be a symmetric (n + 1) x (n + 1) matrix with entries Aij < i, j < n. Define 

n n 

QiA) = AooJ2Au-Y.Al. 

1 = 1 4=1 

Thus Q is a quadratic function on the vector space of symmetric (n + 1) x (n + 1) 
matrices. 

Lemma 1 1. If A>Q then Q{A) > and if A>Q then Q{A) > 0. 

2. If A, B are matrices with Q{A) ~ Q{B) > and if the entries Aqq,Bqq 
are positive then for each s e [0, 1] 

Q{sA + (1 - s)B) > Q{A) , Q{A - B) < 0. 

Moreover, if A ^ B then strict inequality holds. 

To see the first item, observe that we can change basis in R" C R"+^ to reduce 
to the case when the block Aij,l < i, j < n is diagonal, with entries bi say. 
Then if A > we have Aoo&i > Aq^ and so 

with strict inequlality if A > 0. 

For the second item, we just have to observe that Q is induced from a a 
quadratic form of Lorentzian signature on R"+^ by the linear map 

n 

n:A^ {AoQ,^Au,Ao,). 

i=l 

The hypotheses imply that tt{A) and tt{B) are in the same component of a 
hyperboloid defined by this Lorentzian form and the statements follow immedi- 
ately from elementary geometry of Lorentz space. 

Using this Lemma we can deduce the uniqueness of the solution to our 
Dirichlet problem, in any dimension. 

Proposition 1 If (f>o, Cz Ti. then there is at most one solution $ of the equa- 
tion Q{^) ^ € on X X [0, 1] with 1 - Ax^ > for all t and with $(a;,0) = 
0o(a;),*(a;,l) =01 (a;). 
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We show that the functional given by (10) is convex with respect to the 

obvious hnear structure. Thus we consider a 1-parameter family $s = $ + stpj 
with the fixed end points. We have 

ds Jo Jx 

Integrating by parts (just as in the derivation of the geodesic equation) we 
obtain 

^S(<i>.)=^'^(g(<i>,)-e) VdM- 

Suppose that $Oi'i'i are two different solutions, so when s = 0, 1 the term 
g($s) ~ e in the above expression vanishes pointwise. Item (2) in the lemma 
above implies that for s g (0, 1) we have q{^s) — e > 0, with strict inequality 
somewhere. This means that E{^i) > E{^q). Interchanging the roles of $0, $1 
we obtain the reverse inequality, and hence a contradiction. 

One can also prove this uniqueness using the maximum principle. Note 
too that the uniqueness is what one would expect, formally, from the negative 
curvature of the space Ti and the convexity of the functional V . 

4.2 Comparison with the free-boundary hterature 

The author is not at all competent to make this comparison properly. Suffice 
it to say, first, that the problem we are considering is very close to those which 
have been studied extensively in the applied literature. For example, in the 
0-formulation, the condition of prescribing the pull-back of the flux on the free 
boundary is the same as that in the classical problem of the "porous dam" ([1] 
Chapter 8, [8] Chapter 4.4), but with the difference that in that case p is constant 
and there are additional boundary conditions on other boundary components. 
Second, the constructions we have introduced in Section 3 above all appear in 
this literature. The transformation from 6* to taking the harmonic function 
as a new independent variable is called in ([3], Chapter 5) the "isothermal 
migration method" . The transformation from the formulation in terms of 9 to 
that in terms of U is known as the Baiocchi transformation [T], [5], The 
transformation of the free boundary problem for a linear equation to a nonlinear 
Dirichlct problem is used in [10] to derive fundamental regularity results. 

An important feature of the {/-formulation is that it admits a variational 
description. Recall that we are given a function L = max((/)o — 0) on X xR 
and we seek a function U with U > L satisfying the equation A^U = po on 
the set where U > L. This can be formulated as follows. We fix a large positive 
M and consider the functional 

£m{U) = J ^l^xUf + e\d,Uf - poU dpdz, 

over the space of functions satisfying the constraint U > L, where the integral 
is taken over X x [—M, M] in X xR (which, a posteriori, should contain the set 



14 



^Ho,Hi on which U > L). Then the solution minimises £m over all functions 
U > L. This can be used to give another proof of the uniqueness of the solution 
to our problem. It seems likely that it could also be made the basis of an 
existence proof, following standard techniques in the free boundary literature. 
Now recall that our <i>-formulation was based on a variational principle, with 
Lagrangian (10). To relate the two, we consider any function <I> on X x [0, 1] 
with dtdt^ > and define U by the recipe of Section 3. We suppose that 
-M < dt^{x, 0) and dt^{x, 1) < M for all x gX. Then we have 

Proposition 2 The functional £m{U) is 

r I m"^ r f 

E{^)+M / {l-AxM{h-(l>i) + ^\^{^i~M\^dfi+{—+eM) / dfi-e / <j>idfi 

Thus if we fix M and the end points 4>q,4>i the two functionals differ by a 
constant. The central step in the proof is the fact that the integrals 

/ / 52$|Vx$P dfi dt 
Jq Jx 

f f Ax$(at*)2 dn dt 
Jo Jx 

are equal modulo boundary terms. We leave the full calculation as an exercise 
for the reader. 



4.3 The degenerate case 

So far, in this section, we have discussed the case when e > 0. In that case 
the equations we are studying are elliptic. The degenerate case, when e = 0, is 
much more delicate. In fact Chen's main concern in [2] was to obtain results 
about this case, taking the limit as e tends to 0. Chen shows that the Dirichlet 
problem for $, with e = 0, has aC^'^ solution but the question of smoothness is 
open. The formulation of the problem in terms of the function U has particular 
advantages here, because the problem is set-up as a family of elliptic problems, 
and the issue becomes one of smooth dependence on parameters. (This is related 
to another approach, involving families of holomorphic maps, discussed in |12| . 
[7].) We can express the central question as follows. Suppose we have a smooth 
function A on a compact Riemannian manifold X and fix a smooth positive 
function p. Let J be the functional 

J{u) = j -|Vup— /9U. 

For each z G R we set — max(A, z) and minimise the functional J over the 
set of functions u > A^. Suppose we know that there is a minimiser which 
is smooth on the open set ^1^ d X where > A^. Let = {(x, z) : x G l^^} C 
X X R. 

Question 2 In this situation, does vary smoothly with z in Q? 
The interesting case here seems to be when z is a critical value of g. 
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5 Relation with Nahm's equations 

We recall that Nahm's equations are a system of ODE for three functions 
Ti, T2, T3 taking values in a fixed Lie algebra: 

5 = K,Tfc], (19) 

where k run over cyclic permutations of 1, 2, 3. To simplify notation, let us 
fix on the Lie algebra u{n). It is equivalent (at least in the finite-dimensional 
case) to introduce a fourth function Tq and consider the equations 

^ + [T„T,] = [T„T,], (20) 

with the action of the "gauge group" of U'(n)-valued functions u{t): 

Ti uTiU^^ , To 1-^ uTou~^ 

at 

which preserved solutions to (20). (That is, using the gauge group we can 
transform Tq to 0.) The equations imply that 

^(T2 + iTs) = [To + m, T2 + iTs], (21) 
at 

so T2 + iTs moves in a single adjoint orbit in the Lie algebra of GL{n,C). 
Conversely if we fix some B in this complex Lie algebra, introduce a function 
g{t) taking values in GL{n, C) and define skew-Hermitian matrices Ti{t) by 

To + iTi = ^g-\ 
as 

T2+iT^ = 9Bg-\ 

then two of the three Nahm equations are satisfied identically. The remaining 
equation can be expressed in terms of the function h{t) = g*{t)g{t), taking 
values in the space Ti. of positive definite Hcrmitian matrices, which we can also 
regard as the quotient space GL{n, C)/U{n). This equation for h{t) is a second 
order ODE which is the Euler-Langrange equation for the Lagrangian 



m = j \'^^?n + VB{h)dt. 

Here | \u denotes the standard Riemannian metric on Ti,. The function V onH 
is 

VB{h)=TY{hBh-^B*). 
If g is any element of GL{n, C) with g* = h then 



Veih) = \gBg- 



-1|2 
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so Vb is determined by the norm of matrices in the adjoint orbit of B. (See 
[4] for details of the manipulations involved in all the above.) The result in 
[3], mentioned in the introduction to this article, is that for any two points 
ho, hi G Ti. there is a unique solution h{t) to the Euler-Lagrange equations for 
t e [0, 1] with h{0) = ho, h{l) = hi. 

These constructions go over immediately to the case when U (n) is replaced 
by any compact Lie group and GL{n, C) by the complexified group. We want 
to extend them to the situation where U{n) is replaced by the group Q of 
Hamiltionian diffeomorphisms of a surface S with a fixed area form (or more 
precisely, the extension of this group given by a choice of Hamiltionian). The 
essential difficulty is that this group does not have a complexification. However, 
as explained in [6] , [11] , [12] , the space TL of Kahler potentials behaves formally 
like the quotient space Q'^/Q for a fictitious group Q'^. Thus the problem we 
have formulated in Section 2 can be viewed as an analogue of the desired kind 
provided that our potential function V can be seen as an analogue of Vb in the 
finite-dimensional case. 

If we have a path (pt in Ti with (/)o = and a function /3 : E ^ C we can write 
down a differential equation for a one-parameter family (3t which corresponds, 
formally, to the adjoint action of the complexified group Q'^, with the initial 
condition j3o ~ p. The equation has the shape 



The problem is that this evolution equation will not have solutions, even for 
a short time, in general. But if we suspend for a moment our assumption 
that we are working over a compact Riemann surface and suppose that /3 is a 
holomorphic function then there is a trivial solution P-t = (3. So, formally, the 
functional Vp on Ti is given by the norm of j3 with respect to the meaure 



Even if this integral is divergent, the variation with respect to compactly sup- 
ported variations in (p is well-defined, and this is what corresponds to the gra- 
dient of Vb appearing in the equations of motion. Moreover, we can integrate 
by parts to get another formal representation of a functional with the same 
variation 



Now take the compact Riemann surface E to be a 2-torus, and identify the 
space TL with periodic Kahler potentials on the universal cover C. On this 
cover the identity function (3 is holomorphic, and we see from the above that 
the formal expression 
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is analogous to the function Vb in the finite-dimensional case. Of course the 
integrand is periodic and so the integral will be divergent but we can return to 
the compact surface S and consider the well-defined functional 



which will generate the same equations of motion. So we see that, modulo some 
blurring of the distinction between E and its universal cover, the functional we 
have been considering is indeed analogous to that in the finite-dimensional case. 

Using the transformation from the $ equation to the 9 equation, we obtain 
a relation between Nahm's equations for the Hamiltonian diffeomorphisms of a 
surface and harmonic functions on R^. This can be seen in other ways. Most 
directly, we consider three one-parameter families of functions hi{t) on a surface 
S with an area form which satisfy: 



where { , } is the Poisson bracket. We think of these as a one-parameter family 
of maps : S — > R^, and assume for simplicity that these are embeddings, 
with disjoint images. Then it is a simple exercise to show that the equations 
(22) imply that the images hf{Yj) are the level sets of a harmonic function on a 
domain in R^. From another point of view, the geometric structure defined by 
a solution to the <l> equation is an invariant Kahler metric CI = ojq + idd^ 
on S X 5^ X (0, 1) with volume form 



Since drd/S is an S'^^-invariant holomorphic 2-form, what we have is an 5*^ invari- 
ant hyperkahler structure. Then the relation with harmonic functions appears 
as the Gibbons-Hawking construction for hyperkahler metrics. 

The development above is rather limited, since we have only been able for- 
mulate an analogue of our Nahm's equation problem for a single function /3. 
One can go further, and arrive at other interesting free boundary problems. 
Consider for example the case when the surface S is the 2-sphere with the stan- 
dard area form, and the orientation-reversing map cr : S — > S given by reflection 
in the Xi , X2 plane. Now consider maps /3 : S — > C with /3 = (3 o a which are 
diffeomorphisms on each hemisphere. Then the push-forward of the area form 
on the upper hemisphere defines a 2-form pp on C with support in a topological 
disc C C. (The form will not usually be smooth, but will behave like 

c?^^/^ where d is the distance to the boundary of /9(S).) Clearly the form pp de- 
termines /3 up to the action of the tr-equivariant Hamiltionian diffeomorphisms 
of E. Suppose that ft. is a tr-invariant function on E. We can regard this as 
an element of the Lie algebra of Q'^ and consider its action on /3. This is given 
by Ach where h is thought of as a function on C, vanishing outside /3(E). So 
a reasonable candidate for a model of the quotient of the space of maps (3 by 
the action of is given by the following. We consider 2-forms p supported on 




— = {hj,hk} 



(22) 



9? = dTdj3dTdj3. 
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topological discs in C, with singularities at the boundary of the kind arising 
above, and impose the equivalence relation that po ~ Pi if there is a compactly 
supported harmonic function F on C with IS.F = po — pi. 

Now let 9{xi,X2, z) be a harmonic function on an open set fl C R'^, with 
9{xi,X2,z) = 6{xi,X2, —z). Suppose that Vl is diffeomorphic to S*^ x (0, 1), that 
^ = on the inner boundary component Eq and 6 = 1 on the outer boundary 
component Ei. Suppose also that the projections of Eo,Ei to the {xi,X2) 
plane are diffeomorphisms on each upper hemisphere, mapping to a pair of 
topological disc Dq C Di. Then the flux of S/0 on each boundary component 
pushes forward to define a pair of compactly supported 2-forms po,pi on C. 
These are equivalent in the sense above, since po — pi = ^cF for the function 



where the integral is taken over the intersection of the vertical line through 
(a;i,a;2,0) with 17. Our hypotheses imply that -F > 0, and F is supported on 
the larger disc Di. 

The question we are lead to is the following 

Question 3 Suppose Dq C Di are topological discs in C, that pi are 2-forms 
supported on Di and that there is a non-negative function F on C, supported on 
Di, with po — pi — AcF (where the Laplacian is defined in the distributional 
sense). Do pQ, pi arise from a unique harmonic function 9 on a domain in R^, 
by the construction above? 

(For simplicity we have not specified precisely what singularities should be al- 
lowed in the forms pf. this specification should be a part of the question.) 

Hitchin showed in 9J that Nahm's equations form an integrable system. The 
root of this is the invariance of the conjugacy class given by (21), together with 
the family of similar statements that arise from the SO{3) action on the set-up. 
In this vein, we can write down infinitely many conserved quantities for the 
solutions of our equation (11) on the Riemannian manifold Ti.. Let f\ be an 
eigenfunction of the Laplacian Ax, with eigenvalue A > 0. Then we have 

Proposition 3 For any e > 0, if 4't satisfies (11) then the quantity 



does not vary with t. 

This becomes rather transparent in the 6'-formulation, using the fact that 
the function 



satisfies A^K\ = 0. 

The author is grateful to Professors Colin Atkinson, Xiuxiong Chen, Darryl 
Holm and John Ockcndon for helpful discussions. 
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